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A  polymer  electrolyte  membrane  fuel  cell  (PEM  FC)  system  as  a  power  source  used  in  mobile  applications 
should  be  able  to  produce  electric  power  continuously  and  dynamically  to  meet  the  demand  of  the  driver 
by  consuming  the  fuel,  hydrogen.  The  hydrogen  stored  in  the  tank  is  supplied  to  the  anode  of  the  stack 
by  a  fuel  delivery  system  (FDS)  that  is  comprised  of  supply  and  recirculation  lines  controlled  by  different 
actuators.  Design  of  such  a  system  and  its  operation  should  take  into  account  several  aspects,  particularly 
efficient  fuel  usage  and  safe  operation  of  the  stack. 

The  exiting  unconsumed  hydrogen  is  circulated  and  reused  to  increase  the  efficiency  and  at  the  same 
time  maintain  the  humidity  in  the  anode  side  of  the  stack,  thereby  preventing  drying  and  flooding  in 
the  channel  which  can  affect  the  stack  performance.  A  high  pressure  difference  across  a  cell  between  the 
anode  and  cathode  could  cause  damage  on  thin  layers  of  the  cell  components  and  water  imbalance  in  the 
membranes. 

In  this  paper,  we  analyze  a  hybrid  fuel  delivery  system  that  consists  of  two  supply  and  two  recirculation 
lines.  The  major  components  were  a  pressure  regulator,  a  flow  control  valve,  an  ejector,  and  a  blower.  These 
models  were  developed  and  connected  in  order  to  analyze  dynamic  behavior  of  the  fuel  delivery  system. 
Based  on  the  models,  two  control  strategies,  a  decentralized  classic  proportional  and  integral  control  and 
a  state  feed-back  control  were  designed  and  optimized  to  keep  a  constant  pressure  in  the  anode  flow 
channel  and  a  constant  ratio  of  mass  flow  rates  from  recirculation  to  supply  lines.  The  integrated  system 
with  the  two  different  controllers  was  simulated  to  evaluate  its  tracking  and  rejection  performance  at 
different  references  and  disturbances. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  important  technologies  that  must  be  developed 
to  meet  future  power  generation  needs  [1  ].  Applications  in  fuel  cell 
vehicles  have  the  potential  to  substantially  to  reduce  emissions  and 
increase  engine  efficiency  [2].  Among  the  various  fuel  cell  technolo¬ 
gies,  the  polymer  electrolyte  membrane  (PEM)  fuel  cell  will  play  a 
major  role  in  the  future  hydrogen  economy  [3]  and  is  considered 
the  best  candidate  to  replace  the  combustion  engine  because  of  its 
capability  of  high  power  densities,  low  operating  temperatures,  and 
short  start-up  time  [4]. 

Hydrogen  stored  in  a  tank  is  highly  pressurized  in  order  to 
increase  the  volume  density  of  the  fuel.  By  contrast,  the  pressure 
of  supplying  hydrogen  in  the  stack  is  relatively  low;  therefore,  a 
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pressure  drop  is  necessary  to  allow  the  supply  of  fuel  to  quickly 
respond  to  power  demands.  In  addition,  it  is  desirable  for  uncon¬ 
sumed  hydrogen  to  be  diverted  to  a  supply  line  to  increase  the 
efficiency  of  fuel  usage  and  performance  of  water  management. 
These  requirements  are  fulfilled  by  the  fuel  delivery  system  (FDS). 
Basic  functions  of  the  FDS  are  to  regulate  the  flow  rate  of  hydrogen 
supplied  to  the  stack  and  purge  both  inert  gases,  migrated  nitrogen 
from  the  cathode  and  liquid  water  residing  in  the  stack.  Therefore, 
a  supply  line  and  a  purge  line  are  necessary  in  a  fuel  delivery  sys¬ 
tem.  The  supply  line  allows  the  hydrogen  from  the  storage  to  flow 
into  the  anode  of  the  fuel  cell  stack,  where  the  pressure  and  mass 
flow  rate  are  regulated.  The  purge  line  periodically  discharges  the 
hydrogen  from  the  anode  and  at  the  same  time  removes  the  liquid 
water  and  inert  gases  accumulated  in  the  channel. 

A  high  excessive  ratio  of  hydrogen  is  generally  preferred,  which 
reduces  the  response  time  at  an  increased  power  demand,  improves 
efficiency,  and  eases  water  removal;  the  latter  is  particularly 
advantageous  in  water  management  [5,6].  In  these  studies,  the 
unconsumed  fuel  exiting  at  the  outlet  of  the  anode  flow  channel 
was  diverted  to  a  supply  line  by  an  extra  ejector  or  blower  [5]. 
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Nomenclature 

aw  water  activity 

A  area(m2) 

cp  specific  heat  at  constant  pressure(J  kg-1  K-1 ) 

Cd  discharge  coefficient 

d  diameter  (m) 

F  faraday  constant  (96487  C  mol-1 ) 

/  current  (A) 

J  rotational  inertia  (kg  m2 ) 

m  mass  (kg) 

M  molecular  weight  (kg  mol-1 ) 

Ma  Mach  number 

Nceu  number  of  fuel  cell 

p  pressure  (Pa) 

R  gas  constant  (J  kg-1  K-1 ) 

t  time  (s)  or  thickness  (m) 

T  temperature  (K) 

U  speed  (ms-1) 

V  volume  (m3) 

W  mass  flow  rate  (kg  s-1 ) 

y  mass  fraction 

Greek  letters 

oinet  net  water  transfer  coefficient 

y  specific  heat  ratio 

rj  efficiency 

A  water  content 

p  density  (kg  m-3) 

r  torque  (Nm) 

0  relative  Humidity 

0  scaled  mass  flow  rate 

W  dimensionless  heat  parameter 

co  angular  velocity  (rad  s-1 ) 

Subscription 

an  anode 

bl  blower 

bm  blower  motor 

ca  cathode 

dry  dry  state 

ej  ejector 

em  ejector  manifold 

fc  fuel  cell 

fcv  flow  control  valve 

g  mixed  gas  of  hydrogen  and  vapor 

H2  hydrogen 

H20  water 

in  inlet 

l  liquid  water 

Ipr  low  pressure  regulator 

m  membrane 

max  maximum 

memb  membrane 

min  minimum 

mix  mixing  section  of  ejector 

out  outlet 

p  primary  flow  of  ejector 

py  primary  flow  at  secondary  flow  choking  section 
reacted  reacted  hydrogen 
ref  reference 

s  secondary  flow  of  ejector 

sat  saturation 


set 

setting 

sm 

supply  manifold 

St 

stack 

sy 

secondary  flow  at  chocking  section 

V 

vapor 

* 

steady  state 

When  the  load  varies  dynamically,  as  is  usually  found  in  mobile 
applications,  the  FDS  should  supply  the  fuel  quickly,  keeping  vari¬ 
ation  of  the  pressure  in  the  anode  flow  channel  of  the  stack  to  a 
minimum  in  order  to  prevent  a  potential  high  pressure  difference 
across  the  individual  cells  and  subsequent  adverse  effects  on  water 
transport  in  the  membranes  [7]. 

The  dynamics  of  the  FDS  were  studied  by  several  authors  using 
models  of  components  and  systems.  Pukrushpan  [8]  proposed  a 
model  for  a  system  where  a  recirculation  line  and  the  anode  side 
of  the  stack  were  not  considered  in  the  study.  Bao  et  al.  pre¬ 
sented  models  for  a  fuel  supply  line  and  a  recirculation  line  with 
an  ejector  [6],  but  the  recirculation  flow  rate  was  not  controlled 
independently  because  the  ejector  was  in  the  critical  flow  condi¬ 
tion.  Karnik  and  Sun  [7,9],  proposed  a  control-oriented  model  that 
considered  one  supply  line  and  one  recirculation  line  with  an  ejec¬ 
tor  to  circulate  the  exiting  unconsumed  hydrogen.  In  their  study, 
the  recirculation  flow  was  only  controlled  by  the  back  pressure 
valve. 

We  analyzed  a  new  hybrid  fuel  delivery  system  that  increases 
efficiency  of  fuel  usage,  dynamically  supplies  the  fuel  to  the  stack, 
and  at  the  same  time  ensures  safe  operation  of  the  stack.  The 
FDS  consists  of  two  supply  and  two  recirculation  lines.  The  supply 
lines  are  operated  based  on  the  load  demand.  At  a  relatively  low 
load  demand,  the  supply  line  with  a  low  pressure  regulator  mainly 
accounts  for  the  supply  of  fuel.  By  contrast,  the  other  line  with  a 
flow  control  valve  is  used  to  supply  additional  flow  at  a  high  load 
demand.  An  ejector  and  a  blower  serve  to  mix  the  exiting  uncon¬ 
sumed  fuel  with  the  supplying  flow  through  two  recirculation  lines. 
The  ejector  is  a  passive  device,  while  the  blower  is  actively  used  to 
control  the  recirculation  flow  rate.  The  control  objectives  were  to 
maintain  the  hydrogen  stoichiometric  ratio  (SR)  and  the  pressure 
difference  between  anode  and  cathode  using  two  actuators,  a  flow 
control  valve  and  a  blower.  The  schematic  diagram  of  the  hybrid 
FDS  is  shown  in  Fig.  1. 

The  following  sections  are  divided  into  three  parts:  description 
of  the  models  of  major  components  and  the  hybrid  system,  steady 
state  analysis  of  the  FDS,  and  design  of  controls. 

2.  Models  for  components  and  hybrid  system 

The  hybrid  FDS  is  composed  of  manifolds,  an  anode  flow  chan¬ 
nel,  an  ejector,  a  blower,  a  pressure  regulator,  and  a  flow  control 
valve.  Our  modeling  (1 )  assumed  no  contaminants  in  the  fuel  sup¬ 
plied  from  hydrogen  tank  and  no  pressure  drops  along  the  pipes, 
(2)  neglected  spatial  variations,  (3)  used  the  ideal  gas  law  for  all 
volumes,  (4)  assumed  isothermal  conditions  for  manifolds  and  the 
anode  gas  flow  channel,  and  (5)  assumed  that  the  condensed  liq¬ 
uid  water  in  the  control  volumes  is  only  removed  by  the  purging 
operation. 

2.1  Manifolds 

Pipes  and  connections  in  the  FDS  were  approximated  by  two 
manifolds,  an  ejector  manifold  and  a  supply  manifold.  The  ejector 
manifold  represents  the  pipe  connecting  the  flow  control  valve  to 


J.  He  et  al.  /  Journal  of  Power  Sources  185  (2008)  973-984 


975 


Valve 

Fig.  1.  A  hybrid  fuel  delivery  system. 


the  ejector,  and  the  supply  manifold  represents  the  pipe  connecting 
the  ejector  and  blower  to  the  stack. 

Because  of  the  assumption  that  the  gas  in  the  ejector  mani¬ 
fold  is  pure  hydrogen  without  contaminants,  its  pressure  dynamics 
were  described  by  applying  the  mass  conservation  law  in  a  control 
volume; 


saturated  in  the  volume,  the  relationship  between  the  water  activ¬ 
ity  and  humidity  is  given  by; 

0(.)  =  min(aWi(.),  1)  (7) 

By  substituting  Eq.  (4)  and  Eq.  (6)  into  Eq.  (3),  the  dynamics  of 
the  water  activity  are  derived  as; 


dpem 

~dT 


Rh2 T em 

Vem 


(Wfcv 


we 


ej,PJ 


(1) 


da™ 


dt 


Ru2oTSm 

Psat{TSm)VSl 


-(w„ 


!  —  Wv,sm,out) 


(8) 


The  inflow  of  supply  manifold  includes  the  recirculation  gases 
from  the  ejector  and  the  blower,  which  are  a  mixed  state  of  hydro¬ 
gen  and  water  vapor.  Using  the  assumption  that  no  liquid  water 
enters  or  leaves  the  supply  manifold,  the  dynamics  of  the  hydrogen 
and  water  in  the  supply  manifold  are  described  by; 


The  species  mass  flow  rates  at  the  inlet  and  outlet  of  supply  man¬ 
ifold  in  Eq.  (2)  and  Eq.  (8)  are  determined  by  the  ejector  model  and 
anode  inlet  model.  The  inlet  of  the  anode  flow  channel  is  regarded 
as  a  nozzle,  and  the  mass  flow  rate  is  calculated  by  the  nozzle 
equations  [11]; 


dpH2,sm 

dt 


Rh2  Tsm 


(Wh  2,sm,i'n  —  Wh2  ,sm,out) 


(2) 


dm 


H20,  sm 


dt 


—  W v,sm,in  ~  W v,sm,out 


0) 


When  the  mass  of  water  in  a  volume  exceeds  its  saturation  limit, 
liquid  water  will  condenses  in  it.  The  mass  of  saturation  water  vapor 
in  a  volume  is  determined  by; 


mv,sat,(-) 


Psat(,T(.))V(.) 

Kh2oT(.) 


(4) 


where  (•)  represents  the  index  of  volumes  with  vapor,  sm  and  an, 
respectively,  referring  to  the  supply  manifold  and  anode  flow  chan¬ 
nel;  mv,sat  is  the  saturation  mass  of  water  vapor  in  the  volume. 
Psat  is  the  saturation  pressure  (Pa),  given  in  the  function  of  volume 
temperature  (K)  obtained  by  curve  fitting  [10]; 

logio(Psat(T))  =  1.44  X  10“7r3  -  9.18  X  10“5T2 

+2.95  x  1(T2T- 2.18  (5) 


For  a  volume  of  hydrogen  gas  with  water  vapor  and  liquid  water, 
the  water  activity,  aw  and  humidity,  0  are  defined  as  [11  ]; 


«w,(0 


mH20,(Q  . 


and 


+  _  Pv,( .)  . 

* H  Psat^)' 


(6) 


According  to  the  ideal  gas  law,  the  water  activity  is  equal  to  the 
humidity  for  sub-saturated  vapor  in  the  volume.  When  the  vapor  is 


^nozzle  ~  Pu®{z)  (9) 


where 


where  pu  is  the  upstream  pressure,  z  the  pressure  ratio  of 
upstream  and  downstream,  Co  the  discharge  coefficient  of  the  noz¬ 
zle,  and  At  is  the  throat  area  of  the  nozzle;  R  and  y  are  the  gas 
constant  and  specific  heat  ratio  of  the  upstream  gas. 

When  applying  Eq.  (9)  and  Eq.  (10)  to  the  inlet  of  the  anode 
flow  channel,  the  pressure  in  the  supply  manifold  is  the  upstream 
pressure,  the  pressure  in  the  anode  flow  channel  is  the  downstream 
pressure,  and  the  inlet  of  the  anode  flow  channel  is  considered  the 
throat  of  a  nozzle.  The  average  gas  constant  and  specific  heat  ratio 
in  the  supply  manifold  used  in  Eq.  (10)  are  defined  as; 


yg,()  =  KH2yH2,(.)  +  kh2o(i  -yH2,(.)) 

(11) 

Rg,()  =  RH2yH2,(.)  +  rh2o(i  -yH2,(.)) 

(12) 

where 

_  Ph2,()mh2 

H2,()  Ph2,(.)Mh2  +  Psat(T(.))(/)(.  )Mh2o 

(13) 

where  yn2  is  the  mass  fraction  of  hydrogen  in  the  volume,  yg  the 
average  specific  heat  ratio  of  humidified  hydrogen  in  the  volume 
based  on  mass  fraction,  and  Rg  is  the  average  gas  constant. 
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2.2.  Anode  gas  flow  channel 


The  anode  gas  flow  channel  was  considered  an  isothermal 
lumped  volume  with  no  crossover  and  leakage.  Similar  to  the  man¬ 
ifold  models,  the  pressure  and  the  water  activity  in  the  anode  flow 
channel  are  described  as  follows; 


dpu2,an  Rh  -Jan 


dt 


dayman  RH20Tan 


dt 


w  (Wh2  ,an,in  W H2,an,out  ^H2,  reacted)  (14) 

van 

r\Tnn 

(W y^anjn  ~  Wv, an, out  ~  ^v,m)  (15) 


Psat(Tan)Vai 


where  PH2,an  is  the  partial  pressure  of  hydrogen  in  the  anode  flow 
channel.  The  species  mass  flow  rates  entering  the  channel  are 
expressed  as; 


where  subscription  (•)  represents  an  and  ca,  respectively  refer¬ 
ring  to  the  anode  and  cathode. 

Hence,  the  transport  rate  of  vapor  across  the  membrane  from 
the  anode  the  cathode  was  determined  at  a  given  stack  current  and 
water  activities  of  the  flow  channels  at  anode  and  cathode  sides. 

One  control  objective  for  the  FDS  is  to  maintain  a  specified 
small  pressure  difference  between  anode  and  cathode  and  pre¬ 
vent  a  potential  rupture  of  the  membrane.  The  total  pressure  in  the 
anode  flow  channel  is  calculated  by  the  partial  pressures  of  species 
as; 

Pan  =  Ph2  ,an  +  Psat(Tan)(fi  an  (24) 

where  0an  =  min  (l,aw,an)  is  the  relative  humidity  of  gas  in  the 
anode  flow  channel. 


Wh: 2,an,in  ~  WanjnyH2,smi  and  W v,an,in  ~  ^an,m(  1  yH2,sm)  (16) 

where  Wan  in  is  the  mass  flow  rate  entering  anode  flow  channel 
derived  from  Eq.  (9)  and  Eq.  (10),  andyH2,Sm  is  the  hydrogen  mass 
fraction  in  the  supply  manifold. 

WH 2, an, out  and  WVtan,out  in  Eq.  (14)  and  (15)  are  the  mass  flow 
rates  of  hydrogen  and  vapor  leaving  the  anode  flow  channel,  which 
are  derived  from  the  blower  and  ejector  models  in  the  following 
sections.  WH2  jeacted  in  Eq.  (14)  is  the  rate  of  hydrogen  consumed  in 
the  electrochemical  reaction  given  by  [11  ]; 

^reacted  =  Ncell^f^  (17) 

where  NceU  is  number  of  cells  in  a  stack,  Ist  the  stack  current  and  F 
is  Faraday’s  constant. 

Wv,m  in  Eq.  (15)  is  the  transport  rate  of  water  vapor  through  the 
membrane  from  the  anode  to  the  cathode  resulting  from  back  dif¬ 
fusions,  electro-osmotic  drag  forces,  and  pressure  gradients  along 
the  membrane.  The  total  transport  rate  of  water  vapor  is  given  by 
[12]; 

14/  iw  IstM H20 

Wv,m  —  Ot  net N cell - p -  (18) 

where  anet  is  the  net  water  transfer  coefficient  per  proton. 

Because  of  a  small  pressure  gradient  across  the  membrane,  the 
water  transport  by  pressure  gradients  is  negligible,  and  the  net 
water  transfer  coefficient  is  given  by  [11,12]; 

anet  =  nd-  (Xca  ~  Kn)  (19) 

1st  tmlVlmdry 

where  AfC  (cm2)  is  the  active  area,  pm>dry  and  Mmdry  are  the  density 
and  the  weight  per  mole  of  the  dry  membrane,  tm  is  the  membrane 
thickness,  nd  the  electro-osmotic  drag  coefficient,  and  Dw  is  the 
diffusion  coefficient  [11  ]; 

nd  =  0.0029Aj|n  +  0.05A.a„  -  3.4  x  10“19  (20) 

D.-D, exp  (2416(jlj-i))  (21) 

where  Dx  is  a  diffusion  coefficient  related  to  anode  water  content 
defined  by  [12]; 


2.3.  Ejector 


Ejectors  have  been  widely  applied  in  refrigerators  [13,14].  The 
same  techniques  are  used  in  the  recirculation  systems  of  the  FDS, 
which  allow  the  low  pressure  gas  exiting  from  the  anode  to  mix 
with  the  high  pressure  gas  in  the  supply  line  [7,15,16]. 

Based  on  the  one-dimensional  constant  pressure  mixing  theory 
and  the  critical  mode,  the  primary  and  secondary  inlet  flows  at  an 
ejector  are  likely  to  be  choked  [13].  By  Eqs.  (9)  and  (10),  the  primary 
and  secondary  inlet  mass  flow  rates  are  [7]; 


Wej,P  - 


Wej,s  = 


PemAt  ej  /  yw2 


2  x  (n^+DAn^-i) 


y/r^  V  Rh2  V>+  +  1 


PanAsy  /  Yg  ,an 

\fEan  V  R8’an  V  Yg’an  +  1 


( Yg,an+l)/(Yg  ,an~l ) 


(25) 

(26) 


where  r]p  and  ps  are  the  efficiency  for  the  primary  and  secondary 
flow.  At  ej  is  the  nozzle  throat  area  of  the  ejector,  Yg,an  and  Rgtan  are 
the  average  specific  heat  ratio  and  gas  constant  of  the  anode  flow 
channel,  respectively,  and  is  the  hypothetical  throat  area  equal 
to  the  secondary  flow  section  area  where  the  flow  is  choked  [9], 
which  was  determined  by; 

A  _  J  Amix  -Apy,Apy  <  Amfx 

Asy  -  1  n  A  >  A  ■  ^z/) 

{  u’  npy  -  n mix 

where  Amix  is  the  mixing  tube  section  area  of  the  ejector  and  Apy  is 
the  primary  flow  section  area  determined  by  [13]; 

(kh2+i)/[2(kh2-i)1 


(28) 


PexpMapy  \  Yh2  +  7 


1  +  ^ —Mapy 


Ma2  =  + 


(29) 


o 

1 

o 

for 

Xa  <  2 

10-,0(l+2(+n-2)) 

for 

2  <  Xa  <  3 

10_1°(3  —  1.67(A.an  —  3)) 

for 

3  <  Xa  <  4.5 

1.25  x  10“10 

for 

Aa  >  4.5 

(22) 


In  the  above  equations,  Aca  and  Xan  are  the  water  contents  of 
the  membranes  at  the  cathode  and  anode  sides  obtained  from  the 
water  activity  aw  as  [11  ]; 

{  0.043  +  17.81aW  (.)  -  39.85cF  ( }  +36.0cF  ( }  for  aWj(.)  <  1 

14  +  1.4(aW)(.)  -  1)  for  1  <  aW)(  )  <  3  (23) 

16.8  for  aWj(.)  >  3 


where  rjeXp  is  the  coefficient  that  accounts  for  loss  of  primary  flow 
affected  at  the  boundary  [13],  Ma  is  the  Mach  number  of  primary 
flow  at  the  section  where  the  secondary  flow  is  choked.  In  addition, 
the  hydrogen  and  water  vapor  mass  flow  rates  at  the  outlet  of  the 
ejector  are  obtained  by  applying  the  mass  conservation  principle; 

2,ej,out  =  ^ej,p  +  Wej  syu2,an  (50) 

W v,ej,out  —  ^ej,s(l  —  yH2,an)  (51) 

where yn2, an  is  the  hydrogen  mass  fraction  in  the  anode  flow  chan¬ 
nel  calculated  by  the  Eq.  (11 ). 
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2.4.  Blower 

A  blower  used  in  the  recirculation  loop  allows  the  humidified 
hydrogen  to  circulate  from  the  anode  to  the  supply  line.  The  blower 
model  contains  two  parts,  a  blower  map  and  an  electric  motor.  The 
blower  map  is  used  for  the  static  model,  which  determines  mass 
flow  rate  and  blower  efficiency,  while  the  rotational  inertia  of  the 
electric  motor  was  considered  for  the  dynamics.  The  blower  speed 
is  derived  from  the  model. 

The  model  for  the  blower  assumes  that  the  inlet  gas  of  the  blower 
is  the  same  as  that  from  the  return  manifold  and  the  blower  outlet 
pressure  is  equal  to  that  of  the  supply  manifold. 

Because  the  mass  flow  rate  of  the  blower  depends  on  the  inlet 
gas  conditions,  the  corrected  mass  flow  rate  and  angle  velocity  are 
defined  as  [17]; 


Wbc  =  Wbl 


\/Trm/Tref 

Prm/Pref 


and  cobc  = 


coM 


yj Trm/Tref 


(32) 


where  Wbl  is  the  mass  flow  rate  of  the  blower,  Tref  the  reference  tem¬ 
perature  at  288  (I<),  pref  the  reference  pressure  of  one  atmosphere 
pressure,  and  oobi  is  the  blower  angle  velocity.  The  dimensionless 
head  parameter  ^bl  is  given  in  [17]; 


Cp,anTan((Psm/PaniYg'an-')/rg’an  ~  1 ) 

*“ - wm, - 


(33) 


where 


Cp,an  —  cp,H2yH2,an  +  cp,v( 1  ~yH2,an) 


(34) 


Ubl  = 


dbl^bc 

2 


(35) 


where  cPAn  is  the  average  specific  heat  in  the  return  manifold  with  a 
constant  pressure  of  the  humidified  hydrogen,  cPiV  and  cp?h2  are  the 
specific  heats  with  a  constant  pressure  of  water  vapor  and  hydro¬ 
gen,  Ubi  is  the  blower  rotor  tip  speed,  and  db{  is  the  diameter  of 
blower  rotator.  The  scaled  blower  flow  rate  0bi  is  defined  by  [17]; 


&bi 


Wbc 

Pan'll  ^d^Ubi 


(36) 


where 


Pan  — 


Pan 

Rg,anTan 


(37) 


where  pan  is  the  gas  density  in  the  return  manifold.  Then,  the  scaled 
blower  flow  rate  0bi  and  blower  efficiency  rjbi  are  expressed  as 
follows  [17]; 


®bi 


Qi  0bi  +  a2 
^bl  ~  a3 


where  a,  =  axl  +  aaMa,  i  =  1,2,3 


(38) 


rjbl  =  +  b2<Pbl  +  b3  where  b,  =  b,t  +  baMa,  i  =  1,2,3 

(39) 


Ma  = 


Ubi 


,rmRg,rmTrm 


(40) 


where  a  and  b  are  function  parameters,  and  Ma  is  the  Mach  num¬ 
ber  of  the  blade  tip  velocity.  The  parameters  are  derived  by  curve 
fitting  of  the  experimental  data.  The  calculation  results  are  shown 
in  Table  1. 

Based  on  the  static  blower  map  above,  the  mass  flow  rate 
through  the  blower  is  given  as  a  function  of  the  pressures  and  tem- 


Table  1 

Blower  map  function  parameters. 


a 

Value 

b 

Value 

an 

-1.598  x  10“3 

bn 

-3.20  x  104 

au 

2.663  x  nr2 

bn 

1.51  x  105 

<321 

-3.06  x  10-2 

b2\ 

1.78  x  102 

0-22 

-0.174 

b22 

-1.77  x  102 

CI31 

14.6 

£>31 

-2.83  x  10“2 

a32 

-15.7 

b32 

-0.326 

peratures  of  inflow  and  outflow  and  blower  angle  velocity,  which 
is  derived  from  the  blower  motor  as  [8]; 


dT-jbl{rbm  Xb>) 

r  _  CpjrrJrm  (  IPjm  \  Vg,rm 

Tb,~  COb,pb,  \\Prm) 


-l)  Ww 


(41) 

(42) 


Tbm  —  Pbm  n  fabl  kv^bl) 
^ bm 


(43) 


where  Jbi  is  the  rotational  inertia  of  rotator,  r\bm ,  kt,  kv  and  Rbm  are 
motor  constants,  and  ubt  is  the  control  voltage  of  the  blower  motor. 


2.5.  Pressure  regulator 


The  low  pressure  regulator  can  adjust  its  outlet  pressure  to  a 
setting  pressure  by  manipulating  the  mass  flow  rate  through  an 
embedded  nozzle.  The  nozzle  throat  area  is  varied  by  the  outlet 
pressure.  In  this  study,  a  static  pressure  regulator  was  used  due  to 
short  response  time.  The  scaled  mass  flow  rate  <2>/pr  and  the  pressure 
drop  0 ipr  are  defined  as  follows; 


^ Ipr  ~ 


Wipr 
W ipr,  max 


and  >/v  =  P;pr’set  P'Pr’0Ut 

Pref 


(44) 


where  W/pr  is  the  mass  flow  rate  through  the  low  pressure  regula¬ 
tor,  Wiprjnax  the  maximum  flow  rate,  Pipr>set  the  setting  pressure  of 
the  regulator,  Pipr<out  the  outlet  pressure,  and  pref  is  the  reference 
pressure  used  for  scaling,  which  is  equal  to  the  setting  pressure  in 
our  model. 

By  using  the  curve  fitting  method,  the  scaled  mass  flow  rate  at 
the  press  regulator  is  approximated  with  a  polynomial  as; 


^ Ipr  —  C\03  +  C202  +  +  0minjpr 


(45) 


where  <2>min,/Pr  =  W iPr,minlWiPr,max ,  and  Wipr>min  is  the  minimum  con¬ 
trollable  mass  flow  rate.  The  constants  are  obtained  by  a  curve 
fitting  method  using  the  experimental  data  as  shown  in  Table  2. 

Because  the  pressure  and  temperature  of  hydrogen  tank  is 
regarded  as  constant,  the  mass  flow  rate  through  the  regulator  is 
only  a  function  of  the  supply  manifold  pressure  (outlet  pressure). 
In  addition,  it  should  be  noted  that  the  scaled  mass  flow  rate  <2>/pr 
calculated  from  Eq.  (44)  becomes  saturated  at  1. 


2.6.  Flow  control  valve 


The  flow  control  valve  is  regarded  as  a  nozzle  with  a  variable 
throat  area,  and  its  steady  state  behavior  is  approximated  by  a  linear 


Table  2 

Constants  for  the  equation  used  for  the  low  pressure  regulator. 
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function  of  control  input  signal; 

Wfcv  =  UfcvWfcv  ,max  (46) 

where  is  the  control  input  signal  of  the  valve  that  varies  from 
0  to  1,  and  Wfcv  max  is  the  mass  flow  rate  at  the  completely  opened 
throat. 


3.  System  analysis 

The  operation  conditions  of  the  hybrid  FDS  were  divided  into 
two  modes:  low  and  high  loads,  each  dependent  upon  the  output 
power  of  the  stack. 

In  the  low  load  mode,  the  low  pressure  regulator  line  completely 
supplies  the  hydrogen  consumed  by  the  fuel  cell  stack,  while  the 
flow  control  valve  for  the  supply  line  is  closed  and  the  ejector  recir¬ 
culation  line  is  out  of  operation  because  no  fuel  is  being  fed  to  the 
primary  inlet  of  the  ejector.  The  blower  is  the  only  pump  to  cir¬ 
culate  the  excess  hydrogen  exiting  from  the  anode  flow  channel. 
Thus,  the  operation  in  the  low  load  was  simplified  with  a  pressure 
regulator  and  a  blower  recirculation  system.  The  pressure  in  the 
anode  flow  channel  of  the  stack  is  maintained  by  the  low  pressure 
regulator  and  the  blower  allows  adjustment  of  the  hydrogen  recir¬ 
culation  flow  rate.  However,  the  anode  pressure  and  excess  flow 
could  not  be  simultaneously  controlled  by  one  actuator  because  of 
the  coupled  effects  between  the  supply  line  and  the  recirculation 
loop. 

In  the  high  load  mode,  the  low  pressure  regulator  alone  is  not 
able  to  meet  the  fuel  requirements  for  a  high  output  power  of  the 
stack.  The  flow  control  valve  for  the  supply  line  begins  comple¬ 
menting  the  additional  fuel  needed.  Both  the  ejector  and  the  blower 
are  in  operation  to  circulate  the  excess  hydrogen  exiting  from  the 
anode  flow  channel,  where  the  pressure  in  the  anode  flow  channel 
and  hydrogen  flow  rate  in  recirculation  are  controlled  by  two  actu¬ 
ators,  the  flow  control  valve  and  blower.  Because  of  the  complexity 
of  the  high  load  mode  and  its  importance  in  dynamic  responses, 
the  following  investigations  are  focused  on  analysis  of  the  FDS  in 
the  high  load  mode  and  design  of  associated  controllers. 


3.1.  Control  problem  formulation 

The  dynamic  system  model  of  the  FDS  was  obtained  by  con¬ 
necting  the  component  models  given  in  Eq.  (l)-(46).  The  model  of 
the  FDS  includes  six  state  variables:  the  gas  pressure  in  the  ejector 
manifold,  the  hydrogen  pressures  in  the  supply  manifold  and  anode 
flow  channel,  water  activities  in  the  supply  manifold  and  anode 
flow  channel,  and  finally  the  blower  angular  velocity.  The  flow  con¬ 
trol  valve  and  the  blower  are  the  two  actuators  used  for  controlling 
the  FDS.  The  stack  current  is  considered  as  a  disturbance  input  to 
the  plant.  The  resulting  state  equations  for  the  FDS  are; 

r  x  =f(x,  u ,  w) 

X  =  ( Pem »  PH2,srm  aw,sm>  PH2,am  aw,an ,  &>b/)  ( AH\ 

\  ,  xT  \H/ J 

u  =  (ufcv,  UM) 

.  W  =  Jst 


where  x  is  the  state  vector,  u  the  input  vector,  and  w  is  the  distur¬ 
bance  input. 

The  control  objectives  of  the  FDS  are  to  maintain  the  pressure 
difference  permissible  between  the  anode  and  cathode  and  con¬ 
trol  the  hydrogen  flow  rate  for  the  stack  at  dynamically  varying 
loads.  The  hydrogen  SR,  defined  as  a  ratio  of  the  hydrogen  flow 
rate  supplied  to  that  reacted  in  the  anode  flow  channel  was  used 
as  a  performance  variable.  Hence,  the  system  output  is  the  anode 


Table  3 

Model  parameters. 


Symbol 

Value 

Symbol 

Value 

Prc 

1  x  106  Pa 

Q  ,an,in 

0.95 

Trc 

293  K 

NCell 

381 

W fcv.max 

2  x  10“3  kg  s_1 

Van 

5  x  10“3  m3 

Vem 

4  x  10-3  m3 

Tan 

353  K 

Tem 

293  K 

Tst 

353  K 

A t,ej 

8.04  x  10“6  m2 

Afc 

0.028  m2 

Am,ej 

4.07  x  10“5  m2 

Pm 

2  x  103  rrr3 

Op 

0.88 

tm 

5  x 10“4  m 

Os 

0.80 

Mm 

1.1  kg  mol-1 

Oexp 

0.70 

Jbl 

2.6  x  10“3  kgm2 

Pset,lpr 

1.5  x  105  Pa 

kt 

0.15  Nm  A-1 

Wipj-max 

1.753  xlO-3  kgs”1 

K 

0.15  Vs  rad-1 

VSm 

4  x  10“3  m3 

Rbm 

0.82  ohm 

Tsm 

353  K 

tlbm 

0.9 

Aan.in 

7.5  x  10“5  m2 

dbl 

0.15  m 

pressure,  and  the  SR  is  expressed  as  follows; 


fy=g(x,  u,w) 

\y  =  [Pan,  SR]T 


(48) 


where  y  is  the  output  vector,  x,  u  and  w  are  defined  in  the  Eq.  (46), 
Pan  is  the  gas  pressure  in  anode  flow  channel  defined  in  Eq.  (16), 
and  SR  =  Wpi2 ^an,in/^H2 , reacted • 

Because  the  state  equations  above  are  nonlinear,  operating 
points  were  obtained  from  steady  state  analysis,  and  then  the  equa¬ 
tions  were  analyzed. 


3.2.  Analysis  at  steady  state 


According  to  the  law  of  mass  conservation  at  a  steady  state,  the 
mass  flow  rate  of  hydrogen  consumed  by  the  stack  should  be  equal 
to  the  sum  of  the  two  supply  lines  as; 

WH2, reacted =  Wtpr +  Wfcv  (49) 

where  W ^  reacted  is  the  hydrogen  consumption  rate  at  a  steady 
state  and  as  a  function  of  the  stack  current  1st.  W*pr  and  W*cy  are 
the  mass  flow  rates  of  the  low  pressure  regulator  and  flow  control 
valve,  respectively. 

Because  there  is  no  external  humidification  on  the  anode  side, 
the  membranes  of  the  stack  are  considered  a  water  source  or  sink 
for  the  FDS  system.  At  a  state  of  equilibrium,  the  water  flow  rate 
cross  the  membrane  should  be  zero,  according  to  the  mass  con¬ 
servation  principle.  When  it  is  assumed  that  the  membrane  at  the 
cathode  side  is  completely  humidified  in  the  high  load  mode,  the 
water  content  at  cathode  side  becomes  the  maximum  value  of  16.8. 
The  water  transport  across  the  membrane  by  back  diffusion  is  bal¬ 
anced  by  the  electro-osmotic  drag  force.  As  a  result,  the  net  water 
transfer  coefficient  becomes  zero  at  steady  state.  Thus,  the  corre¬ 
sponding  equilibrium  water  activity  in  the  anode  flow  channel  at 
a  steady  state  for  a  given  stack  current  is  calculated  by  setting  the 
net  water  transfer  coefficient  to  zero  in  Eq.  (19). 

A  resulting  typical  Ist-aWf(m  curve  at  the  water  equilibrium  state 
is  depicted  in  Fig.  2,  and  the  membrane  parameters  used  in  the 
calculation  are  listed  in  Table  3. 

As  shown  in  Fig.  2,  two  different  values  of  the  equilibrium  water 
activities  in  anode  flow  channel  at  a  given  stack  current  existed; 
which  value  results  depend  upon  the  previous  state.  For  example, 
when  the  stack  current  is  increased  from  0  to  400  A,  water  activity 
decreases  along  the  line  AB,  BC  and  finally  the  line  after  C.  By  con¬ 
trast,  when  the  stack  current  is  decreased  from  400  A  to  0,  the  water 
activity  increases  along  the  line  CDA.  The  liquid  water  may  be  con¬ 
densed  in  anode  flow  channel  at  the  low  current  range  of  less  than 
about  12  A  (dry-out  current),  where  the  equilibrium  water  activity 
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Fig.  2.  Water  activity  versus  stack  current  at  water  equilibrium. 


is  larger  than  1.  When  the  stack  current  is  larger  than  the  dry-out 
current,  the  water  activity  becomes  less  than  1,  and  no  liquid  water 
exists  in  the  anode  flow  channel.  In  the  high  load  mode,  the  rela¬ 
tionship  of  water  activity  and  stack  current  at  steady  state  can  be 
approximated  using  a  linear  polynomial. 

At  a  given  stack  current,  values  for  the  state  variables  and  control 
variables  can  be  obtained  by  steady  state  analysis.  First,  the  anode 
water  activity  is  calculated  by  the  water  balance  analysis  shown  in 
Fig.  2,  and  the  humidity  and  hydrogen  mass  fraction  in  the  anode 
flow  channel  can  be  obtained  from  Eq.  (6)  and  Eq.  ( 13 )  at  the  desired 
pressure  in  anode  flow  channel  given. 

Then,  the  water  vapor  mass  flow  rate  leaving  the  anode  flow 
channel  is  calculated  using  the  following  relationship; 


w*  _  -  w*  _ 

vvv,an,out  ~  H2,an,out  y 


1  —  V*  1  —  V* 

_ fcp*  1  U/l/*  U2,an 

-  (SR  1  )WH2, reacted  " 


H2,an 


y\ 


H  2, an 


(50) 


Because  water  activity  is  at  a  steady  state,  the  water  vapor  mass 
flow  rate  entering  the  anode  flow  channel  should  be  equal  to  that 
leaving  the  channel.  The  mass  flow  rates  entering  and  leaving  the 
anode  flow  channel  are; 

Wan,in  =  SR*WH2, reacted  +  WZ,an,out’ 

and  W*„j0Ut  =  (SR*  -  1  )W*2  reacted  +  w;„  (51 ) 


Likewise,  the  hydrogen  mass  fraction  in  the  supply  manifold  is 
obtained  by  the  mass  flow  rates  entering  the  anode  flow  channel 
given  by; 


SR*W; 


^H2,sm 


H2 ,  reacted 


w* 


(52) 


In  addition,  the  supply  manifold  pressure  is  calculated  by  the 
inverted  nozzle  equations,  Eqs.  (9)  and  (10),  from  the  mass  flow  rate 
entering  the  anode  channel  and  hydrogen  mass  fraction.  The  mass 
flow  rate  of  the  low  pressure  regulator  is  attained  by  the  supply 
manifold  pressure.  Both  the  mass  flow  rate,  W^,  and  the  control 
input  signal  of  the  flow  control  valve,  u|cy,  are  derived  by  Eq.  (48) 
and  Eq.  (45). 


At  a  steady  state,  the  mass  flow  rate  at  primary  inlet  of  ejector  is 
equal  to  that  of  flow  control  valve.  The  pressure  of  the  ejector  man¬ 
ifold  is  derived  from  Eq.  (24).  The  mass  flow  rate  at  the  secondary 
inlet  of  the  ejector,  W*.  s  is  obtained  by  the  ejector  model  in  the  part 
2,  which  yields  the  blower  mass  flow  rate; 

w*bl  =  W*„i0Ut  -  w*j  s  (53) 

The  angular  velocity  of  blower  was  obtained  using  the  inversion 
of  static  map  of  blower  in  Section  2  by  the  mass  flow  rate  and  the 
pressure  differences  of  supply  manifold  and  anode  flow  channel. 
Thus,  the  control  input  u*,  is  determined  by  the  steady  state  form 
of  Eq.  (17)  and  Eq.(21). 

Flence,  the  reference  control  variables  ujcv  and  u*t  are  a  function 
of  control  objectives,  pressure  in  the  anode  flow  channel,  hydrogen 
SR,  and  stack  current  from  the  above  analysis.  Because  the  pres¬ 
sure  in  the  cathode  channel  is  a  function  of  the  stack  current,  the 
pressure  in  anode  flow  channel  also  is  expressed  as  a  function  of 
the  stack  current  [6].  For  an  FDS  designed  for  a  100  kW  fuel  cell 
stack,  the  reference  pressure  in  the  anode  flow  channel  is  given 
by; 


p*n  =  -0.04583 Js2t  -  7.917/st  +  150000  (54) 

where  p*n  is  the  reference  value  of  the  pressure  (Pa)  in  the  anode 
channel,  and  Ist  is  the  stack  current  in  the  range  of  300-400  A  for 
the  high  load  mode. 

The  reference  hydrogen  SR*  is  assumed  to  be  1.5  for  different 
stack  currents.  Other  parameters  are  given  in  Table  3. 

The  operating  points  of  the  system  were  found  using  tools  avail¬ 
able  in  Matlab.  Ten  different  stack  currents  were  selected  and  the 
results  for  the  steady  state  analysis  are  shown  in  Fig.  3.  As  seen  in 
Fig.  3(a)  and  (b),  the  control  signals  for  the  flow  control  valve  and 
the  blower  are  nearly  linear  to  the  stack  current. 

The  flow  rate  of  consumed  hydrogen  and  the  mass  flow  rates 
entering  and  leaving  the  anode  flow  channel  increases  with  the 
increasing  stack  current.  Conversely,  the  mass  flow  rate  of  the  low 
pressure  regulator  decreases  because  of  the  increased  pressure  in 
the  supply  manifold.  In  addition,  the  recirculation  mass  flow  rate 
of  the  ejector  secondary  inlet  decreases.  Consequently,  the  flow 
control  valve  for  the  supply  line  becomes  the  primary  fuel  sup¬ 
ply  line  and  the  blower  becomes  the  primary  circulation  pump  for 
hydrogen  recirculation  as  stack  current  increased. 

As  shown  in  Fig.  3(d),  the  supply  manifold  pressure  increases 
and  the  pressure  in  the  anode  flow  channel  decreases,  as  defined 
in  Eq.  (53),  as  the  stack  current  increases.  In  addition,  the  pres¬ 
sure  drop  between  the  inlet  and  outlet  of  the  anode  flow  channel 
increases.  Thus,  as  stack  current  increases,  higher  parasitic  power  is 
required  in  the  blower,  which  was  dissipated  to  increase  the  mass 
flow  rate  in  the  recirculation  loops  and  to  overcome  the  pressure 
drop  along  the  anode  flow  channel. 

The  calculated  six  state  variables  of  the  system  at  the  steady 
state  are  shown  in  Fig.  4.  As  the  stack  current  increases,  the  hydro¬ 
gen  partial  pressure  in  the  supply  manifold  increases  as  shown 
in  Fig.  4(a),  and  that  in  the  anode  flow  channel  decreases.  In 
addition,  water  activity  in  the  anode  flow  channel  and  in  the 
supply  manifold  decreases,  but  activity  in  the  anode  flow  chan¬ 
nel  is  higher  than  the  other,  as  shown  in  Fig.  4(b).  The  ejector 
pressure  and  the  blower  angle  velocity  are  almost  linear  to  the 
magnitude  of  the  stack  current,  as  shown  in  Fig.  4(c)  and  (d).  All 
the  six  state  variables  were  determined  using  a  curving  fitting 
method  based  on  the  stack  current  and  its  relationship  given  in 
Fig.  4(a)-(d). 


980 


J.  He  et  al.  /  Journal  of  Power  Sources  185  (2008)  973-984 


Fig.  3.  (a)  Reference  flow  control  valve  input  at  a  steady  state,  (b)  Reference  blower  control  voltage  input  at  steady  state,  (c)  Mass  flow  rates  at  steady  state,  (d)  Pressures  at 
steady  state. 


3.3.  Static  feed-forward  (SFF) 

In  the  steady  state  analysis,  the  reference  control  variables  are 
expressed  as  a  function  of  the  stack  current  at  a  desired  anode  pres¬ 
sure  as  given  in  Eq.  (30)  and  the  assumed  hydrogen  SR*  =  1.5.  Based 
on  that  analysis,  the  curves  in  Fig.  3(a)  and  (b)  are  approximated  by 
the  following  two  equations  used  to  obtain  control  signals  for  the 
flow  control  valve  and  the  blower; 


U*fcv  =  2.108  X  10“3/st  -  0.1391 

(55) 

u*bl  =  0.931 1st  -  137.9 

(56) 

Based  on  these  equations,  a  static  map  was  created  as  a 
static  feed-forward  (SFF)  block  shown  in  Fig.  5  and  used  for 
determining  the  inputs  for  the  FDS.  The  outputs  are  the  pres¬ 
sure  in  the  anode  flow  channel  and  the  hydrogen  stoichiometric 
ratio. 


4.  Design  of  linear  controller 

The  SFF  had  some  drawbacks  including  unexpected  distur¬ 
bances,  parameter  uncertainties,  and  approximation  errors  in  the 
map  data  that  caused  deviations  in  the  steady  state  response.  For 
example,  map  data  that  represents  the  relationship  between  the 
voltage  for  the  motor  and  the  mass  flow  rate  is  nonlinear.  The  SFF 
control  of  the  blower  is  based  on  a  curve-fitted  equation  that  inter¬ 
polated  the  inverse  map  data.  Any  change  of  parameters  like  aging 
effects  might  cause  deviations  in  the  characteristic  curve  for  the 
blower. 

Linear  feed-back  controls  are  the  first  choice  to  avoid  instability 
of  the  system.  For  design  of  the  controls,  the  nonlinear  FDS  should 
be  linearized  at  an  operating  point,  which  was  selected  at  the  stack 
current,  Jst  =  350  A.  The  resulting  control  signals  for  the  flow  control 
valve  and  the  blower  are  0.59  and  187  V,  respectively.  The  state  form 
of  the  linearized  FDS  is; 

8x  =  A8x  +  Bu8u  +  Bw8w  .  . 

8y  =  C8x  +  Dw8w 
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Stack  Current  lst(A) 


Fig.  4.  (a)  Hydrogen  pressure  at  steady  state,  (b)  Water  activities  at  steady  state,  (c) 
Pressure  of  ejector  manifold  at  steady  state,  (d)  Angle  velocity  of  blower  at  steady 
state. 


follows; 

G(s)  =  kp  +  y-  (58) 

where  kp  and  lq  are  the  proportional  and  the  integral  gains  of  PI 
controllers,  respectively. 

The  gains  are  determined  using  a  design  tool  given  in  Matlab  for 
the  SISO  system.  The  chosen  PI  controller  gains  for  UfCV-pan  system 
are  /<p, pan  =  6.79  x  10-5  and  kIpan  =  5.22  x  10-5  s_1,  while  those  for 
Ubi- SR  system  are  kPt sr  =  4.58  and  7</sr  =  915  s-1 . 


4.2.  State  feed-back  control 

Another  feed-back  control  is  state  feed-back  control  with  inte¬ 
gral,  which  can  stabilize  the  system  by  feed-backing  outputs  and 
state  variables.  The  resulting  full  state  feed-back  control  (SFB)  with 
integral  is  shown  in  Fig.  7,  where  two  control  loops  are  involved. 
The  inner  loop  serves  to  enhance  the  dynamics  of  the  response  by  a 
state  feed-back  controller,  while  the  outer  loop  compensates  for  the 
steady  state  tracking  error  via  an  integral  controller.  In  addition,  the 
outputs  of  the  SFF  are  feed-forwarded  to  the  control  loop  to  reduce 
the  transient  response  time  needed  for  settling  to  a  steady  state 
from  a  previous  state. 

Based  on  the  linear  model  and  the  block  diagram  above,  the 
control  variable,  u,  is  expressed  as  a  function  of  the  feed-forward 
output  u,  the  perturbation  of  the  state  variables  8x,  and  the  output 
of  the  integrator  q; 

u  =  u*  —  K8x  —  KjQ  (59) 

q=y*-y  (60) 

where  I<  and  K\  are  the  state  feed-back  and  integral  gain  matrices, 
respectively,  and  are  optimized  using  the  linear  quadratic  regulator 
(LQR).  The  cost  function  for  the  LQR  is  shown  as; 

poo 

J=  (8yTQy8y +  qTQ/q  +  8urR8u)dt  (61) 

Jo 

The  cost  function  includes  an  output  vector,  8y,  that  is  a 
function  of  the  disturbance  term,  8w,  as  shown  in  Eq.  (56). 
To  simplify  the  cost  function,  a  new  output  vector  is  defined 
as; 

8/  =  C8x  (62) 


Substituting  8y'  to  8y  in  the  Eq.  (60),  the  cost  function  yields  a 
new  form; 

where  8{  )  =  ()-(  )0  refer  to  the  perturbation  from  operating  points, 

( )o  are  the  variables  at  operating  points,  and  A,  Bu ,  Bw,  C  and  Dw  are 
the  system  matrices. 

The  FDS  with  the  controller  is  depicted  in  Fig.  7,  where  the  inte¬ 
grator  q  is  considered  as  a  new  state  variable; 

4.1  Classic  proportional  integral  control  q  =  y  *  -y  =  -8y  =  -C8x  -  Dw8w  (64) 


/ 


(8xCTQyC8x  +  qTQiq  +  8utR8u)  dt 


(63) 


Although  the  outputs,  pressure  in  anode  flow  channel,  and 
hydrogen  SR  are  dependent  on  both  inputs,  the  FDS  can  be  con¬ 
sidered  as  two  independent  single-input  and  single-output  (SISO) 
systems,  UfCV-pan  system  and  u^-SR  system,  to  design  a  diagonal 
controller  as  shown  in  Fig.  6.  The  use  of  a  diagonal  controller  for  a 
system  with  a  multi-input  and  multi-output  (MIMO)  is  referred  to 
as  a  decentralized  controller.  Thus,  the  anode  pressure  is  controlled 
by  the  flow  control  valve,  while  the  hydrogen  SR  is  controlled  by 
the  blower. 

The  decentralized  PI  controller  is  implemented  by  two  inde¬ 
pendent  PI  controllers,  whose  transfer  function  is  given  as 


Neglecting  the  perturbation  of  disturbance  signal,  8w,  the 
expanded  state  equation  integrating  Eqs.  (56)  and  (63)  is  given 
as; 


8x 

Q 


A  0 
-C  0 


8x 

+ 

V 

8  u  =  A 

8x 

+  B8u 

Q 

0 

_Q 

(65) 


The  optimization  gains  of  control  laws  in  Eq.  (58),  which  min¬ 
imize  the  cost  function  in  Eq.  (62)  for  an  infinite  period,  is  given 
as; 

[k  k/]  =  j?-1bp 


(66) 
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Fig.  5.  Static  feed-forward  control. 


Fig.  6.  Decentralized  PI  controller. 


where  P  is  obtained  from  the  solution  of  the  Algebraic  Riccati  equa¬ 
tion; 


pA+Atp  +  q_-  pbr~'btp  =  o 


(67) 


where 


linearization  model.  Selected  Qy,  Qj  and  R  are  as  follows; 


Qy  = 


4  x  10“6 


0.4 


,  Q/  = 


6  x  10 


,-5 


32 


R  = 


16 


4  x  10“ 


(69) 


/  2.595  x  10_1  1.293  xl0“4 

1^-2.591  xlO4  -8.401  x  10"1 


Q  = 


CTdyC 


Qi 


Thus,  the  matrices  of  the  controller  gains  are  computed  by  solv¬ 
ing  Eq.  (65)  and  (66)  as  follows; 


-1.561 


9.758  x  10“ 


1.404  x  lCr  1.373  x  10“ 


-2.435 

3.018 


10“ 


1.865  x  10“ 


(68) 


I<i  = 


-1.021  x  10“2 
-1.234  xl0“3  -1.090 


2.984  -1.803  xlO3 

4.3.  Simulation  results  and  analysis 


(70) 

(71) 


When  the  weighting  matrices  Qy,  Qi  and  R  are  given,  the  control 
gains  I<  and  K\  are  obtained  by  the  LQR  algorithm  based  on  the 


Fig.  8  shows  the  transient  response  of  the  pressure  and  the 
hydrogen  SR  at  a  step  current  change.  The  operating  point  for  the 


Fig.  7.  State  feed-back  control. 
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Fig.  8.  (a)  Step  change  of  stack  current,  (b)  Pressure  change  in  anode  flow  channel, 
(c)  Hydrogen  stoichiometric  ratio  change. 


stack  current  Ist  was  350  A,  and  a  current  step,  AJst  =  10A,  was 
applied  after  1  s  as  shown  in  Fig.  8(a).  Fig.  8(b)  shows  the  pressure 
change  responses  in  the  anode  flow  channel  for  the  three  controls. 
The  settle  time  of  the  SFB  took  only  about  2  s,  while  the  SFF  and 
PI  controls  took  more  than  5  s.  The  SFF  was  not  able  to  reach  the 
desired  steady  state  and  stayed  with  an  error,  while  the  SFB  and  PI 
control  followed  a  reference  pressure  and  finally  reached  the  steady 
state  value. 

Similarly,  the  responses  of  hydrogen  SR  change  are  shown 
in  Fig.  8(c).  The  settle  time  at  disturbance  rejection  of  the  SFB 
was  shorter  than  that  of  the  SFF  and  PI  control.  In  addition, 
the  PI  control  produces  a  higher  overshoot  in  the  SR  response 
than  that  produced  by  the  SFB,  while  the  SFF  control  could 
not  remove  the  steady  state  error.  Consequently,  the  SFB  con¬ 
trol  outperformed  the  other  controls  by  its  highly  dynamic 
response. 

Fig.  9  shows  the  responses  of  the  FDS  at  a  multi-step  stack  cur¬ 
rent  that  varied  from  300  A  to  400  A  with  a  period  of  5  s  as  shown 
in  Fig.  9(a).  Fig.  9(b)  shows  the  flow  control  valve  signal  responses 
under  different  control  strategies.  The  control  signal  of  the  SFB 
shows  a  high  overshoot,  while  it  took  longer  for  the  PI  control  to 
attain  a  stable  state.  In  addition,  the  voltage  for  the  blower  at  the 
SFB  control  also  had  a  high  overshoot,  but  the  settle  time  was  the 
shortest,  as  shown  in  Fig.  9(c)-(e)  show  the  transient  responses  of 
the  pressure  in  the  anode  flow  channel  and  the  hydrogen  SR,  where 
the  SFB  yields  the  best  performance  in  controlling  the  pressure  and 
the  hydrogen  SR.  Fig.  9(f)  shows  the  response  of  the  stack  current  on 
the  water  activity  in  the  anode  flow  channel  where  no  significant 
differences  for  the  three  controls  exist.  The  water  activity  by  all  of 
the  controls  was  smaller  than  1  and  no  liquid  water  was  formed  in 
the  flow  channel.  As  a  result,  no  flooding  occurred  in  the  anode  flow 
channel  at  the  given  multi-step  change  of  stack  current  in  high  load, 
which  makes  justified  the  assumption  of  neglecting  purge  opera¬ 
tion.  Because  the  control  requirement  for  the  pressure  in  the  anode 
flow  channel  is  stricter  than  that  for  the  hydrogen  SR,  the  SFB  is  the 
choice  for  the  FDS. 
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Fig.  9.  (a)  Multi-step  stack  current,  (b)  Flow  control  valve,  (c)  Blower  control  voltage, 
(d)  Anode  pressure,  (e)  Hydrogen  stoichiometric  ratio,  (f)  Anode  water  activity. 


5.  Conclusion 

For  this  study,  we  analyzed  a  hybrid  fuel  delivery  system  for 
a  PEM  fuel  cell  power  system  that  was  designed  to  meet  require¬ 
ments  of  mobile  applications.  The  system  consists  of  two  supply 
and  two  recirculation  lines  with  an  ejector,  a  blower,  and  valves 
that  specifically  promote  efficient,  dynamic  and  reliable  operation 
of  the  fuel  supply.  First,  models  for  individual  components  of  the 
system  were  developed  and  used  to  construct  the  hybrid  system. 
Based  on  the  models  for  the  system,  different  control  strategies 
were  designed  and  their  performances  compared.  The  major  find- 
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ings  of  this  study  are  summarized  as  follows: 

1.  The  conditions  for  the  operation  of  a  hybrid  fuel  delivery  sys¬ 
tem  that  depended  upon  the  magnitude  of  the  load  current  were 
analyzed.  We  found  that  the  system  could  be  operated  by  a  low 
pressure  regulator  and  blower  for  a  low  load.  The  pressure  in  the 
anode  flow  channel  and  the  recirculation  mass  flow  rate  were 
strongly  coupled  and  could  not  be  controlled  simultaneously  by 
one  actuator  in  low  load  mode.  By  contrast,  for  a  high  load,  the 
system  needs  to  be  controlled  to  improve  the  tracking  behav¬ 
ior  of  the  fuel  during  sudden  changes  of  the  stack  current.  From 
steady  state  analysis,  both  the  pressure  in  anode  flow  channel 
and  the  flow  rate  of  the  recirculation  were  operable  by  the  flow 
control  valve  and  blower  for  a  high  load.  Finally,  the  pressure  in 
anode  flow  channel  and  mass  flow  rate  of  recirculation  could  be 
separately  and  independently  controlled  at  a  high  load. 

2.  According  to  the  comparison  of  the  performance  of  different  con¬ 
trol  strategies  that  include  the  SFF,  the  PI,  and  the  SFB,  the  SFB 
control,  based  on  the  linearized  model  of  the  FDS,  demonstrated 
the  best  performance. 

The  major  contributions  of  this  paper  are  ( 1 )  the  development  of 
a  dynamic  model  for  components  of  a  hybrid  fuel  delivery  system, 
(2)  the  design  of  a  decentralized  PI  controller  and  a  multi-variable 
linear  controller,  SFB  controller  that  dynamically  controls  the  recir¬ 
culation  mass  flow  rate  while  variation  of  the  pressure  in  the  anode 
flow  channel  is  suppressed,  and  (3)  the  analyses  of  integrated  sys¬ 
tems  along  with  their  controls. 

In  real  situations,  the  purging  process  is  necessary  to  remove 
inert  gases  that  build  up  in  the  recirculation  loops  as  well  as  the 
liquid  water  that  forms  in  the  anode  flow  channel.  Effects  of  the 
purging  process  on  the  performance  of  FDS  will  be  included  in  a 


future  study.  Measurement  of  all  states  is  not  realistic,  so  reduc¬ 
tion  of  the  states  will  be  carried  out  using  the  observer  for  the  SFB 
control. 
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